Nodeless Superconductivity in the Noncentrosymmetric Superconductor MgioIrigBi^ 
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We measured the resistivity, diamagnetization and low temperature specific heat of the newly discovered 
noncentrosymmetric superconductor Mg\oIriqBi(,. It is found that the superconducting gap has an s-wave sym- 
metry with a value of about Aq ~ 0.6 meV, and the ratio Ao/hgT^ « 1.86 indicates a weak coupling for the 
superconductivity. The correlations among the normal state Sommerfeld constant y„, the slope -dpoHc2(T)/dT 
near and the condensation energy Ec are all consistent with the weak coupling picture. The separated phonon 
contribution from the specific heat shows that the conduction electrons of the Ir atoms interact most strongly 
with the vibrations of themselves, instead of with that of the light element boron. 

PACS numbers: 74.20.Mn,74.20.Rp, 74.25.Bt, 74.70.Dd 
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The study on superconductivity in noncentrosymmetric ma- 
terials has attracted growing efforts in recent years iH 0, S 
3, m. For most superconductors, the atomic lattice has a 
centrosymmetry, therefore the system is inversion symmet- 
ric. The orbital part of the superconducting order parameter 
has a subgroup which is confined by the general group of the 
atomic lattice. Due to the Pauli's exclusion rule and the parity 
conservation, the Cooper pair with orbital even parity should 
have anti-parallel spin state, namely spin singlet, while those 
having orbital odd parity should have parallel spin state, i.e., 
spin triplet. If a system lacks the centrosymmetry, the triplet 
pairing may be instable leading to a mixture of singlet and 
triplet pairing. Theoretically novel features are anticipated in 
the noncentrosymmetric system|6]. A nodal gap structure has 
been observed in Li^PfsB showing the possibility of triplet 
pairing, while due to weaker spin-orbital coupling|0 [st], the 
nodal gap has not been observed in a material Li2Pdj,B with 
similar structure. It is thus highly desired to investigate the 
paring symmetry in more materials with noncentrosymmetric 
structure. 

The newly discovered superconductor MgmIr\<)B\(, (here- 
after abbreviated as MglrB )[9] with superconducting transi- 
tion temperature ~ 5 K is one of the rare materials which 
have the noncentrosymmetry. This material has a space group 
of / - 43ot with large and complex unit cells each has about 45 
atoms. To some extent it resembles the system LiziPt, Pd)T,B 
since it has alkaline metals (Li, Mg), heavy transition ele- 
ments (Pd, Ft, Ir) and the light element boron. Theoretically 
it was shown that the major quasiparticle density of states 
(DOS) derives from the d-orbital of the heavy transition el- 
ements. In this paper we present a detailed investigation and 
analysis on the superconducting properties, such as the energy 
gap, pairing symmetry, electron-phonon coupling strength and 
condensation energy etc. in MglrB. Our results suggest that 
the superconductivity in this system is of BCS type with an 
s-wave gap symmetry and a weak coupling strength. 

The samples were prepared in two steps starting from pure 
elements of Mg (98.5%), Ir (99.95%) and B (99.999%) us- 
ing a standard method of solid state reaction. Appropriate 
mixtures of these starting materials were pressed into pellets, 
wrapped in Ta foil, and sealed in a quartz tube with an at- 



mosphere of 95% Ar/5% H2. The materials were then heated 
at 600 ° C and 900 ° C for 40 min. and 80 min., respec- 
tively. After cooling down to room temperature, the samples 
were reground and mixed with another 20% of Mg, then they 
were pressed into pellets and sealed in a quartz tube with the 
same atmosphere as used in the first step. In this process the 
sample was heated up to 900 ° C directly and maintained for 
80 min. The synthesizing process here is similar to the pre- 
vious work reported by the Princeton group |9] but still with 
some differences. For example, we used Mg powder instead 
of flakes to make the mixture more homogeneous. In addition 
the pressure in the sealed quartz tube may rise to nearly 4 atm 
at 900 ° C, which may considerably reduce the volatilization 
of Mg during the synthesis. The resistivity and the AC sus- 
ceptibility were measured based on an Oxford cryogenic sys- 
tem (Maglab-Exa-12). The specific heat was measured on the 
Quantum Design instrument PPMS with temperature down to 
L8 K and the PPMS based dilution refrigerator (DR) down 
to 150 mK. The temperatures of both systems have been well 
calibrated showing consistency with an eiTor below 2% in the 
temperature range from 1.8 K to 10 K. 

The x-ray diffraction (XRD) patterns taken on these sam- 
ples show a single phase with very small amount of impu- 
rity which is comparable to that reported previously [9]. After 
the first round of synthesizing, the superconducting transition 
inspected by the AC susceptibility occurs at about 5 K with 
a relatively wide transition. However, after the second step, 
the transition moves to about 3.7 K with a sharper transition 
width. This indicates a sensitive dependence of Tc on the Mg 
content. In Fig. 1(a) we show the temperature dependence 
of resistivity under different magnetic fields. It is seen that 
the transition width (1% - 99%p„ ) is about 0.2 K. By ap- 
plying a magnetic field the transition shifts to lower tempera- 
tures quickly with a rather low slope -dfioHc2(T) / dTlr^ ^ 0.3 
T/K. Using the Werthamer-Helfand-Hohenberg relation|fl^ 
jU()//c2 = -Q.69dfii)Hc2(T)/dT\T^,Tc, we get the upper critical 
field fJ.()Hc2 - 0.77 T. The AC susceptibility is shown in Fig.l 
(b) revealing a similar behavior as the resistive transition. It 
is interesting to note that the value of the slope dfii)He2{T) / dT 
found here is lower than that in the earlier report[9] showing 
the tunability of superconducting properties in this system. 
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FIG. 1: (color online) Temperature dependence of resistivity (top) 
and magnetic susceptibility (y" and^') (bottom) under different DC 
magnetic fields. It is clear that the DC magnetic field makes the tran- 
sition shift parallel to low temperatures, manifesting a field induced 
pair breaking effect. The inset in top panel shows the resistive tran- 
sition in a wide temperature regime at ij.qH = OT. 
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FIG. 2: (color online) Raw data of specific heat plotted as C/T vs. 
T^. All filled symbols represent the data taken with the DR based on 
the PPMS at various magnetic fields. The open squares show the data 
taken with the PPMS at zero field. The thick solid line represents the 
normal state specific heat which contains both the phonon jpi, and 
the electronic contributions. 



This will be justified in the following analysis. In any case, it 
is safe to conclude that the normal state Sommerfeld constant 
ranges from about 19.6 to 41.7 mJ/molK^. Further analysis 
suggests that the real y„ is close to the upper bound of the 
experimental values, i.e., y„ 41. 7 mJ/molK-^. 

Next we can have an estimation on y„ . In MglrB, the elec- 
tronic conduction is dominated by the 5d band electrons of 
Ir atoms. The DOS at Ef given by the LDA band structure 
calculation[12] is about N{Ef) - 5.51/eVspin. Assuming a 
electron-phonon coupling constant in the system, one 
has 



Shown in Fig. 2 are the raw data of the specific heat. The 
open squares represent the data taken with PPMS, while all 
filled symbols show that taken with the DR. Both sets of data 
coincide very well at zero field. With increasing the magnetic 
field the specific heat jump due to the superconducting tran- 
sition moves quickly to lower temperatures leaving a back- 
ground which is consistent with that above Tc at zero field. 
This provides a reliable way to extract the normal state spe- 
cific heat as shown by the thick solid line since the normal 
state can be described by C/T - y„ + f5T^, where the first and 
the second terms correspond to the normal state electronic and 
phonon contribution, respectively. From the data it is found 
that y8 = 2.03 mJ/molK'^ and y„ =41.7 mJ/molK'. In low tem- 
perature region (~ 0.2 K) the superconducting state exhibits, 
however a residual value yo ~ 22.1 mJ/molK^ indicating a 
non-superconducting fraction of about 53%. This high value 
of non-superconducting fraction is however difficult to be re- 
garded as due to an impurity phase with completely different 
structure as MglrB since the XRD data is quite clean. We thus 
suggest that the superconductivity depends sensitively on the 
relative compositions among the three elements and some re- 
gions without superconductivity have the chemical composi- 
tion and even the structure close to the superconducting phase. 
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Using the band structure value of N{Ff), we have y,, = 
25.98 (1 -I- Ae-ph){mJ I molK^). Taking the lower bound of 
the experimental value j„ - 19.6 mJ/molK^ implies an un- 
physical value - - 0.25. Taking however, the upper 

bound of the experimental value y,, = 41.7, we get A = 0.6. 
Therefore it seems that the real value of y„ is close to the up- 
per bound of the experimental values. We can also use an 
alternative way to estimate y^. In the dirty limit, for a type-II 
superconductor, one has ! 1 ill 
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whereA = 3.81e/jT^kB=0.0081{T/K)(mD.cm)-HmolK^/mJ) 
for MglrB. Using -dfiQHa{T)ldT\T^, ^ 0.3T/K, p„ = 0.816 
inQcm, and taking rj - 1 for the weak coupling case, we have 
y„ = 45 mJ/molK^ which is also close to the upper bound of 
the experimental value. 

In the raw data shown in Fig. 2, one can see an upturn of 
y - C/T in the very low temperature region. This upturn 
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TABLE I: Fitting parameters of Schottky anomaly. 
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is known as the Schottky anomaly, induced by Hfting the de- 
generacy of the states of the paramagnetic spins. We tried a 
two level (S=l/2) model to fit the low temperature data but 
found a poor fitting together with an extremely large Lande 
factor g in the Zeeman energy g^sHefj, where p.B is the Bohr 

magneton, H^ff - ^H- + //^ is the effective magnetic field 
which evolves into Hgff - Hq, the crystal field at zero exter- 
nal field. In MglrB the most possible paramagnetic centers 
are from l/"^ (S-5/2) or Ir^^ (S=2). The system energy due 
to Zeeman splitting in a magnetic field is nl3ll 

Esch = ^ Eiexpi-EilksT)! ^ exp(-Ei/kBT), (3) 

where £, = Mjg/^BHeff and Mj = -S, -Sh-1, S-1, S. 
The specific heat due to the Schottky eff'ect is thus Csch = 
(n/kB)dEsch/dT, where n represents the concentration of the 
paramagnetic centers. For S-5/2 and S - 2 the calculated 
results are very close to each other, therefore we show only 
the fit with 5=5/2 correspoding to /r'*^ (six levels). This 
method allows us to deal with the data at zero and finite fields 
simultaneously. It is known that the Schottky term should be 
zero at T = K. In the superconducting state, the total specific 
heat can be written as: Cmt = Cnons + Ce + Cph + Csch with 
Cnons - JqT the contribution of the non-superconducting re- 
gions, Ce is the electronic part. In the zero temperature limit 
only the contribution of the non-superconducting part is left. 
Applying a magnetic field gives rise to a finite value Ay^, to 
Ce = JeT due to the presence of vortices. Practically, in order 
to fit the Schottky term, we first remove the phonon contribu- 
tion Cph - /3T^, then vertically move the experimental data 
downward with a magnitude yo = 22. 1 mJ/molK^ and a field 
induced vortex term AjeiH) as shown in Table-I. Four sets 
of data after this treatment and the corresponding fits to the 
Schottky term at the fields of /jqH = 0.0, 0.1, 0.2, 0.6 T are 
shown in Fig. 3. It is clear that low temperature upturn can be 
well described by the Schottky effect. The results yielded by 
the fitting are summarized in Table-I. One can see that when 
field is beyond 0.5 T which is close to the upper critical field, 
we take the total normal state value y„ as the removed back- 
ground which leads to a perfect fitting to the Schottky term as 
shown in Fig. 3(d). 

After successfully removing the Schottky term and the con- 
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FIG. 3: (color online) The calculated Schottky anomaly (red solid 
line) compared with the electronic specific heat (symbols) in the low 
temperature regime, where the phonon term CpiJT, field induced 
term Ay^ and the nonsuperconducting term C„o„ JT have been re- 
moved. The detailed values used for the Schottky anomaly and the 
electronic contributions are given in Table-I. 



tribution from the non-superconducting region, we get the 
pure contribution from the superconducting regions (as shown 
in Fig. 4). Note that here we used the value 41.7 mJ/molK^ as 
the normal state Sommerfeld constant y„ if it would contain 
only the pure superconducting phase. One can see that the low 
temperature part is flattened out below about 0.8 K when the 
field is zero. Actually this flattening is already visible in the 
data shown in Fig. 3(a) before the Schottky term is removed. 
Furthermore it can also be justified by the requirement of en- 
tropy conservation. Since the Schottky term gives only very 
small contribution in the high temperature region (above 1 .5 
K here), if je had a power law instead of a flat temperature 
dependence here, the entropy would be clearly not conserved 
yielding a large negative entropy. This is of course unreason- 
able. In Fig.4 we present together the theoretical curve for je 
calculated using the weak coupUng BCS formula 
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where ^ = ^J£^ + A^iT). In obtaining the theoretical fit 
we take the implicit relation A(T) derived from the weak cou- 
pling BCS theory for an s-wave superconductor and use the 
gap Ao and Tc as two trying parameters. The theoretical curve 
fits the experimental data very well leading to an isotropic gap 
value Ao = 0.6meV and = 3.1 5K. The ratio Ao /kBT^ = 1.86 
obtained here is quite close to the prediction for the weak cou- 
pling limit (Ao/^bTc = 1.76). This is self-consistent with the 
conclusion derived from the estimation on y„. In addition, the 
specific heat anomaly at Tc is ACe/ynTlr,. ~ 1-54 being very 
close to the theoretical value 1.43 predicted for the case of 
weak coupling. The inset in Fig.4 shows a field induced part 
Aje- In an s-wave superconductor Aje is mainly contributed 
by the vortex cores and a linear relation Aye °^ Hjn/HdiO) is 
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s-wave fit: T_, = 3.75 K 
A„ = 0.6 meV 
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FIG. 4: (color online) All symbols show the temperature depen- 
dence of the electronic specific heat in the superconducting state 
with the contributions from phonon, Schottky anomaly and the non- 
superconducting fraction removed. The blue solid line shows the 
BCS fitting curve at zero field. The inset shows the field dependence 
of the specific heat in the zero temperature limit. Note that in this 
figure the normal state Sommerfeld constant y„ has been scaled up 
to 41JmJ/molK- as marked by the horizontal line. 




FIG. 5: The condensation energy calculated from the specific heat. 
Inset (a) shows the entropy in the normal and superconducting state. 
Plotted in inset (b) is the difference of the entropy between the normal 
and superconducting state. 
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For a BCS s-wave superconductor, a - I. Taking y„ = 
Al.lmJImolK^ and Aq - 0.6 meV, we found a value of Ec ~ 
154 mJ/mol which is very close to the experimental value 142 
mj/mol. 

Now we get down to the electron-phonon coupling. From 
the normal state value we have derived that /3 ~ 2.03 
mJ/molK'^. Using the relation ©/, ^(l2n'^kBNAZ/5f3y^\ 
where Na - 6.02x10^-' the Avogadro constant, Z=45 the num- 
ber of atoms in one unit cell, we get the Debye temperature 
®D{MgIrB) - 350.5 K. It is known that the Debye tempera- 
ture for metallic Ir is about 36 meV (420 K)[12]. While the 
crystalline boron has a very high Debye temperature ®b- 100 
meV (1200 K) il2ll . The rather low Debye temperature found 
in our experiment @o{MgIrB) - 350.5 K is close to that of Ir 
metal, this may suggest that the conduction electrons from the 
5d band of Ir couple most strongly with the vibrations of Ir 
themselves. However, as we stressed before, it seems that the 
Tc can be tuned to higher values by changing the relative com- 
positions among the three elements Mg, Ir and B. This may 
enhance the electron-phonon coupling and/or the quasiparti- 
cle DOS at Ef. The basic parameters and properties derived 
in this work provide a playground for the future study in this 
interesting system. 

In summary, analysis on the low-temperature data in 
MglrB finds a s-wave pairing symmetry with a gap in the 
weak coupling regime. The conduction electrons interact pri- 
marily with the vibrations of the Ir atoms leading to a weak 
coupling strength A ^ 0.6. 
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anticipatedlT?] in the low field region with A^/Ef <kT Tc- 
This linear relation is well demonstrated by the data below 0.4 
T, indicating another evidence of s-wave pairing symmetry. 

In the following we try to estimate the superconducting 
condensation energy Ec- In calculating Ec we get the entropy 
difference between the normal state and the superconducting 
State hy Sn - Ss - J,j (jn - 'Ye)dT', then Ec is calculated 

through Ec - Jj. (S„ - S s)dT . The data of S„ and Ss as 
well as the diff'erence between them are shown in inset (a) and 
(b) of Fig. 5, respectively. The main frame of Fig. 5 shows the 
temperature dependence of the condensation energy Ec which 
is about 142 mJ/mol. This value can actually be assessed by 
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